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Abstract We derive an analytical expression that re- 
lates the breadth of a streaked photoelectron spectrum 
to the group-delay dispersion of an isolated attosecond 
pulse. Based on this analytical expression, we introduce 
a simple, efficient and robust procedure to instantly ex- 
tract the attosecond pulse's chirp from the streaking 
measurement. We show that our method is robust against 
experimental artifacts. 

Key words attosecond - streaking - electron - tra- 
jectory 



1 Introduction 

The characterization of isolated attosecond pulses has 
played an important role in the development of attosec- 
ond science [2[2 ,3 . The generation and application of 
ever shorter attosecond extreme- ultraviolet (XUV) pulses 



[H[S] relies on knowledge of their time-domain properties, 
which can be obtained by means of attosecond streaking 
measurements [MZ]- 

So far, the main functions of attosecond streaking are 
(i) to characterize the field of an attosecond pulse and (ii) 
to temporally resolve a physical process on the attosec- 
ond scale. Here, we are concerned with the former appli- 
cation of attosecond streaking, that of characterizing an 
attosecond pulse. Much effort has been exerted on the 
development of methods for extracting physical informa- 
tion from the streaking measurement [8,9,10 , with the 
current state-of-the-art being the FROG retrieval algo- 
rithm [II] . The FROG algorithm has already been used 
to characterize the shortest attosecond pulses [5] , and to 
uncover a measured delay of 20 as between photoemis- 
sions from the 2s and 2p sub-shells of neon It is 
relatively robust [T3] , and provides a wealth of informa- 
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tion about the temporal characteristics of the attosecond 
and laser fields [14]. 

However, the application of FROG to attosecond streak- 
ing requires quite stringent experimental requirements, 
such as a sufficient amount of recorded spectra, with 
a delay step between them on the order of the attosec- 
ond pulse's duration. These experimental parameters be- 
come unwieldy as the duration of attosecond pulses ap- 
proaches the atomic unit of time. Moreover, the FROG 
algorithm is a somewhat complicated numerical opti- 
mization procedure, whose output (the attosecond field 
and the laser field) is not transparently related to the 
input (the set of streaked spectra). Thus, errors in the 
reconstructed pulses are difficult to interpret due to the 
FROG algorithm's black-box nature. Although FROG 
provides a complete characterization of the attosecond 
XUV field, the duration of the attosecond pulse is the 
primordial quantity that will be interrogated as attosec- 
ond streaking continues to expand beyond its original 
scope into various research fields. 

In this article, we introduce a simple and robust method 
for quantifying the chirp of an attosecond pulse based on 
an analytical formula we derive from laser-dressed pho- 
toelectron trajectories. Using this formula, we develop 
a method that directly evaluates the attosecond pulse's 
group-delay dispersion from a sequence of streaked spec- 
tra, which in turn sets the pulse's duration provided 
its spectrum is known. Our method avoids the strin- 
gent experimental conditions required for the attosec- 
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ond FROG technique, and provides accurate results with 
very few electron spectra in a matter of seconds. We be- 
gin this article with the derivation of the analytical ex- 
pression for the change in photoelectron bandwidth due 
to the streaking effect, and then introduce our method 
with a numerical example. All quantities are expressed 
in atomic units unless otherwise stated. 

2 Classical electron trajectory analysis of the 
streaking effect 

Let us first consider an attosecond XUV pulse with elec- 
tric field Fx(t) given by 

F x (t) = |F x (*)|e i(fix ' +¥,x(t)) , (la) 
= \Pit 2 + ^fot 3 + ..., (lb) 

where the spectrum of the attosecond pulse is centered at 
fix with small variations in frequency due to the higher- 
order temporal phase ipx(t). The attosecond pulse launches 
electron trajectories that are parameterized with an ini- 
tial time t as well as an electron energy e — p 2 /2. Due 
to the attosecond pulse's finite bandwidth, we consider 
the energy e as an independent variable, while the in- 
dependent variable t is a result of the finite duration 
of the attosecond pulse. Thus, the set of trajectories is 
described by a time-energy distribution with respect to 
{t,e}. 

The final energy eg of an electron, launched at some 
moment t in a continuum permeated by a near-infrared 
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(NIR) laser field, is then 



£S = 



i^e + wxW)-^*)) 



(2a) 



s-pA^{t)+ l -Al{t)+(l-^-).: :s U). (2b) 



where we define the instantaneous frequency wxif) = 
ipx.it) due to the chirp of the attosecond pulse, and A^it) 
is the vector potential of the laser field. Since the change 
in frequency over the temporal profile of the attosecond 
pulse is much smaller than the central frequency fix, 
the last term in (|2b[) is comparably small and can be 
dropped, leading to the simple relation es ~ e — pAi,(t) + 
Aj j (t)/2 for the shift of the photoelectron spectrum. 

It is known [7] that the spectral shift alone is not suffi- 
cient to obtain information about the attosecond pulse's 
chirp because the final energy es is hardly sensitive to 
the temporal phase fxit) of the attosecond pulse. The 
main manifestation of the attosecond chirp in the streak- 
ing measurement is the change in breadth of the streaked 
photoelectron spectrum. To describe this effect, we in- 
terpret (|2aj) as a mapping of the initial time and energy 
of an electron trajectory to a final energy (e.g. measured 
at the detector). To describe the effect of chirp, it is use- 
ful to consider small changes des in the final energy with 
respect to small changes in the initial energy de and time 
dt of the trajectory. The total differential of (|2al) is then 



de s « 1 



Mt) 



L M\^ x ( t)+pFL (t)) dt + de), (3) 

where we again neglect the small terms containing ujx(t). 
The temporal phase of the attosecond pulse appears in 



as Px{t) — pxit), which defines the chirp of the 
attosecond pulse. We have also introduced the electric 
field of the laser pulse F^{t) — —Ai,(t). Thus, the chirp 
of the attosecond pulse and the electric field of the laser 
pulse both influence the spread in final energies resulting 
from the streaking effect. 

To proceed further, we interpret the effects of the 
NIR field on the time-energy distribution of electron tra- 
jectories, as described by in a straightforward man- 
ner. Initial inspection of ([3]) shows that the NIR field im- 
parts an additional energy sweep of pF^if) to the photo- 
electron, resulting in a total chirp j3§{t) — j3x(t)+pFi,(t). 
Furthermore, the NIR field re-scales the energy spread 
by a factor (1 — Ai,(t)/p). As a result, both the NIR 
electric field Fi,(t) and the NIR vector potential A^it) 
have a role in modifying the breadth of the photoelectron 
spectrum. 

In order to account for the effects of the streaking 
field, we recall that the attosecond electron wave packet 
can be viewed as a replica [6] of the attosecond pulse 
Fxit). We model this photoelectron replica as 



X(t) = c~5(*/^x) 2 e i(ect+i/3xt 2 ) ; 



(4) 



where ec is the central photoelectron energy. Naturally, 
since the electron trajectories are launched by the at- 
tosecond pulse, the duration tx of the electron wave 
packet [U] should be nearly the same as that of the at- 
tosecond pulse; and as \(t) is a replica of Fxit), its chirp 
fix is the same as that of the attosecond pulse. For sim- 
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plicity, we assume fix to be constant and we also assume 
that the attosecond pulse is shorter than any relevant 
time scale of the NIR field, so that Fl(£) and A^(t) are 
evaluated at the central time to of the attosecond pulse. 

Now, in order to include the effects of the streaking 
field, we first consider the shift of the photoelectron spec- 
trum due to Al (to ) and the change in bandwidth due to 
the chirp induced by -FL(to)- To this end, we modify the 
wave packet's central energy Ec = Pq/2 and chirp j3x as 
follows: 

ec — > £s = ec - PcA L (t ) + -Al(t ) (5a) 
/8x— »• Ps=Px(t)+ Pc F L (t). (5b) 

With these substitutions, the streaked photoelectron wave 
packet is modeled as 

Xs (t)=e^(*/rx) 2 e i(es*+^st 2 )_ ( 6 ) 

To obtain an expression for the bandwidth of the streaked 
photoelectron spectrum, we note that the streaked pho- 
toelectron spectrum is just a Fourier-transform of xs(t) 
[T5] . Since the streaked wave packet is a Gaussian, the 
Fourier transform of xs(t) can be carried out analyti- 
cally, yielding the following expression for the bandwidth 
of the streaked spectrum: 



S S (t ) = 




(7a) 

7s(*o)=7x+(|r) PcFM, (7b) 

where 5x and 7x represent the bandwidth and group- 
delay dispersion (GDD) — defined as the second deriva- 
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tive of the spectral phase — of the attosecond pulse. The 
quantity r]x — \j (^x^) + (^x7x) 2 is the attosecond 
pulse's time-bandwidth product, with a Fourier-limited 
time- bandwidth product ij^ (rj^ =1/2 for a Gaussian 
spectrum). The quantities rx, Sx and <5s(to) are all taken 
as standard deviations of their respective distributions. 
According to (0, 7x determines the width <5s(to) of the 
streaked spectrum as a function of to- Provided that the 
characteristics of the field-free spectrum (Ox, $x and 
r]^) as well as those of the laser field (Ai,(t) and F^(t)) 
are known, 7x remains the only free parameter. 

In writing ([7]) , we also explicitly included the energy 
re-scaling pre-factor (1 — Ai,(to) / pc) ■ Similar but less 
general expressions for the streaked photoelectron band- 
width were previously derived in [6lll3| from the semi- 
classical expression for streaking |16| . These expressions 
consider photoionization at the zero-crossing of the vec- 
tor potential, A^ito) = 0, where there is no spectral 
shift but only a change in spectral bandwidth due to the 
NIR field. These expressions therefore do not contain 
the bandwidth re-scaling factor (1 — A-^(to)/pc), which 
is needed to accurately represent the bandwidth of the 
streaked spectra at arbitrary delay times to, when the 
NIR field simultaneously shifts the photoelectron spec- 
trum and changes its bandwidth. 

Although was deduced assuming a Gaussian wave 
packet, it actually applies to more general pulse shapes 
owing to the fact that the relation r/ x = (jlx^ + ^xTx 
holds for arbitrary spectra with a constant GDD (see 



M = 



(8) 



where the sums range over the XUV-NIR delays tj . The 
goal of our procedure is to find 7x that best reproduces 
the measured curve 5-^ (to) according to model ((?]). 

As an example, we consider the case of a non-Gaussian 
~ 226 as XUV pulse. This pulse has a constant GDD of 
~ 5885 as 2 . However, since its spectrum (Figure [TJ-b) is 
irregular (jiyp ps 0.5515), i.e. it is asymmetric and con- 
tains some fine structure, its chirp /3x is time-dependent. 



The direct evaluation of attosecond chirp from a streaking measurement 

Appendix A). The following section presents numerical 
examples in further support of this claim. g 



3 A method to directly extract the attosecond 
chirp from a set of streaked photoelectron 
spectra 

Equation serves as the basis for our method to ex- 
tract the attosecond chirp from a streaking measure- 
ment. Our procedure is very straightforward: we evaluate 
the first moments (eg) of the streaked spectra to obtain 
the laser field's vector potential A^(t), which in turn 
gives us the laser's electric field Fi,(t). We also compute 
a curve <5g M ^(to) of standard deviations of the measured 
streaked spectra as a function of the XUV-NIR delay 
to- Lastly, we find the attosecond chirp 7x — the only 
free parameter in ([7]) — which minimizes the discrepancy 
between the widths £ g (to) obtained from the set of 
streaked spectra and those given by the model 0. To 
compare these two, we define a figure of merit 
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Fig. 1 Panels (a) and (b) show sets of 101 streaked photo- 
electron spectra evaluated by solving the TDSE, using streak- 
ing fields with <j>o = and <f)o = tt/2, respectively. Panel (c) 
shows the attosecond pulse's spectrum (solid line) and phase 
(dotted line), while panel (d) displays its temporal intensity 
profile (solid line) and temporal phase (dotted line). 

The streaking field is a NIR pulse given by 



A L (t) = A cos 4 (t/T L ) sin(w L t + 4>o) 



(9) 



with tl ~ 5.743 fs, yielding a 3 fs full width at half maxi- 
mum (FWHM) duration, ojl ~ 2.355 rad/fs correspond- 
ing to a central wavelength of 800 nm and with Aq ps 
—0.41915 a.u., giving a peak intensity of 20 TW/cm 2 . For 
this example, we consider carrier-envelope phase values 
of 0o = (Figure QJa) and (j) Q = tt/2 (Figure [TJb) . 

The simulated streaking measurements, shown in Fig- 
ure [TJa and Figure [l]-b, are composed of a sequence 
of streaked spectra computed for different delays be- 
tween the XUV and NIR fields by propagating the time- 
dependent Schrodinger equation (TDSE) using a split- 
step FFT scheme. The Hamiltonian is that of a single 
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electron in one dimension, assuming a soft-core poten- 
tial with an ionization energy W ~ 59 eV. 

The results of our analytical chirp evaluation (ACE) 
procedure, applied to the spectrograms shown in Figures 
Q}a and[T]-b, are shown in Figures [5] and [3] In both cases, 
we have applied ACE to different subsets of streaked 
spectra, by considering a varying number N of spectra 
about the central delay value to = 0. 

For the case 4>o = 0, Figure [2}a shows a false-color 
plot of the figure of merit M as defined in ©. Darker 
areas correspond to a smaller value of M. When too few 
spectra are considered, Figure [2}a shows a local min- 
imum near 7x = 12 500 as 2 which disappears as more 
spectra (N > 13) are considered. Nonetheless, Figure 
[5J-b shows that we recover the exact GDD (the dashed 
line) from the global minimum to within ~ 4% with as 
few as three spectra. As N increases, the global min- 
imum eventually stabilizes around the red dashed line 
representing the exact GDD, and ACE converges nearly 
to the exact value 7x = 5885 as 2 . Figure [2}c shows that 
the model ([7} reproduces the correct curve flg(to) for the 
exact GDD. In contrast, we found that the attosecond 
FROG retrieval [H| fails to converge when fewer than 
25 spectra are included, for which it recovers a GDD 
7x = 5740 as 2 . 

For 0o = tt/2, Figure |3J-ashows that the figure of 
merit has only one minimum as a function of GDD. This 
minimum quickly converges to the correct GDD as more 
spectra are considered in the evaluation, as displayed in 
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Figure |3}b, and is already accurate to within 0.7% for 
N = 13 spectra. Figure 0}c shows that the model 
once again reproduces the correct curve (hollow circles) 
of streaked breadths for the exact GDD 7 = 5885 as 2 . 

The main advantage of the ACE procedure is that it 
requires very few spectra. As long as A^fto) is properly 
sampled by the delay step between the spectra, there is 
enough information for ACE to recover the GDD of the 
attosecond pulse. In contrast, FROG requires the delay 
step to be on the order of the attosecond pulse's dura- 
tion. To illustrate this point, we apply ACE to a subset of 
the spectra shown in Figure [I}a andQ]-b. Specifically, we 
consider 17 spectra over the interval [— 2fs, 1.84 fs] (con- 
taining 1.5 cycles of the streaking field), with a delay 
step of 240 as between them, i.e. a third of the original 
spectra in [— 2fs, 2fs]. Even with so few spectra, ACE 
still recovered accurate GDD's of 6150 as 2 and 6110 as 2 
for <fio — and <fro — tt/2, respectively. On the other 
hand, FROG failed to converge to anything meaningful 
in both cases, most likely because the delay step was too 
large. 

To further demonstrate ACE's robustness against a 
non-Gaussian spectrum, we consider a clipped version of 
the XUV spectrum shown in Figure[T]-c, for which we re- 
move energy components above 175 eV. Experimentally, 
such a sharp edge in the XUV spectrum might result 
from the beam's transmission through a metallic filter. 
Using the clipped XUV spectrum, we compute sets of 
101 streaked photoelectron spectra, with the same pa- 
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number of spectra considered XUV-NIR delay t (fs) 



Fig. 2 The analytical chirp evaluation (ACE) is applied to 
the streaking example shown in Figure [T}a. Panel (a) is a 
false-color logarithmic plot of the figure of merit M, defined 
by ([SJl, versus the number of spectra (N) considered for the 
ACE procedure. Panel (b) plots the retrieved GDD (squares) 
at the global minimum of M as a function of N . In panels (a) 
and (b), the dotted red line represents the exact GDD. Panel 
(c) shows the energy es (dotted line) and breadth 5g (to) 
(solid line) evaluated from the streaked spectra. The hollow 
circles represent the breadths 5s (to) computed from Q with 
the exact 7x = 5885 as 2 . 

rameters as those displayed in Figure [T}a and [l}b. In 
spite of this heavy clipping, ACE recovers GDD's of 
5940 as 2 and 5960 as 2 for (f>o = and <fio — tt/2, respec- 
tively. In comparison, FROG recovers a GDD of 5620 as 2 
for both </> = and </>o = n/2. 

As previously mentioned, these examples assume a 
constant GDD over an irregular spectral distribution, re- 
sulting in a chirp ftx that depends on time. Since expres- 
sion — which is at the core of the ACE procedure — 
assumes a constant chirp in time, then the chirp pa- 
rameter /?x is interpreted as the average chirp over the 
attosecond pulse's temporal profile. Conversely, if a non- 
uniform GDD was considered, then ACE would have re- 
covered the average GDD over the spectral profile. 
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number of spectra considered XUV-NIR delay t (fs) 



Fig. 3 The analytical chirp evaluation (ACE) is applied to 
the streaking example shown in Figure [T]-b. The data shown 
here are presented in the same manner as in Figure [21 
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Fig. 4 Panels (a) and (b) shows sets of streaked spectra, 
computed by adding Poisson noise to those of Figures [TJa 
and [T]-b, respectively. ACE recovers GDD's of 5590 as 2 and 
6040 as 2 from the spectra in panels (a) and (b), respectively. 

As an additional verification of ACE's robustness, we 
investigate the effect of noise in the streaked spectra. To 
this end, we add noise to the sets of 101 spectra shown 
in Figure QJa and QJb. We assume that the number of 
counts n in a spectral bin follows a Poisson distribution 
P(n;ii) = /x Tt e _ ' J /n!, with an expectation value /i pro- 
portional to the spectral intensity (we set fi = 1 for the 
peak of the spectrogram, corresponding to a very low 
count rate). From these considerations, we compute the 
noisy spectra which are shown in Figure |4l 

Even under such nefarious conditions, ACE recovers 
accurate values of the GDD: 5480 as 2 from the spectro- 
gram shown in Figure 0Ja, and 5830 as 2 from the one 
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in Figure HJ-b. In comparison, FROG recovers GDD's of 
6250 as 2 and 6100 as 2 , respectively. This example demon- 
strates that ACE can tolerate very noisy spectra, and 
moreover that it is robust against errors in the vector 
potential A^(to), as determined from the streaked spec- 
tra. 



4 Conclusion 

In conclusion, we have derived a general analytical ex- 
pression ([7]) for the change in spectral breadth due to 
the streaking effect by considering the trajectories of a 
photoclcctron ejected by an isolated attosecond pulse in 
a laser held. We have used this equation as a basis for a 
method to directly extract the attosecond chirp from a 
sequence of streaked spectra. In contrast to the attosec- 
ond FROG retrieval, the ACE procedure does not re- 
quire streaked spectra to be recorded with a delay step 
on the order of the attosecond pulse duration: it only 
requires the delay step to properly sample the streak- 
ing field. This alleviates many of the experimental con- 
straints related to the current approaches to characterize 
isolated attosecond pulses. In addition, the ACE proce- 
dure is simple to implement, robust against experimen- 
tal artifacts, and fast — taking seconds to execute and 
requiring very few (< 10) streaked spectra. This makes 
ACE ideal for real-time diagnostics in attosecond streak- 
ing measurements. 
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A A relation between duration, bandwidth and 
dispersion for arbitrary spectra 

The following is a proof of the general relation 



2x2 



r 2 = r 2 + 7 2 5 



(10) 



between the duration t, the Fourier- limited duration tq, 
the bandwidth 8 and the group-delay dispersion (GDD) 
7 for a pulse with an arbitrary spectrum and a constant 
GDD; r, To and 5 are taken as standard deviations of 
their respective distributions. 

Let us first define spectral and temporal profiles as 



f(u) = / (w)e^ 
1 



/(*) 



(11a) 



f(u)fJ ut du = F- 1 [f(w)](t). (lib) 



We assume, without lack of generality, that f(u)) and 
f(t) are centered around uj = and t = 0, respectively. 

The duration r is defined as the standard deviation 
of fit), which is the square- root of the variance 



t 2 \f(t)\ 2 dt 



\-itf(t)\ 2 dt (12) 



dt. 



In the following derivation, the prime symbol (""') de- 
notes differentiation with respect to the argument and 
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the pulse is normalized according to \f(t)\ 2 dt = Using (fT5"]) and (|TT)) . we may now represent the du- 

J^°oo l/( w )| 2( l w — 1- Assuming 7 is frequency-independent, ration r as 
then (fTTaj) implies f(w) = /g(w)e5T" 2 +i 7 /(w). Insert- 



ing this expression for /'(w) into the rightmost-hand-side 
of (fF2"1) . we obtain 

/oo 
\f(t)\ 2 dt (13a) 
-OO 

OO 

+ 7 / (l*(t;7)/'(*)+n*;7)/'*(*))dt, 
7(t; 7 )=J- 1 [/oMe i7w2 ](*). (13b) 

In analogy to (fF2")l . the bandwidth- limited duration 
is given by 



dt. 



(14) 



Now, I(t;~f) and /^(u^es 7 " are Fourier transforms of 
each other. Thus, from Parseval's theorem, we have 



\I(t ;l )\ 2 dt= / /uMc^ 2 dc (15) 



dw = 7q, 



where (fT4"f in combination with Parseval's theorem was 
used for the last equation on the RHS of (fT5|) . 

The bandwidth (5 is given, also in analogy to (|12D , as 



<5 2 = / c 2 |/H| 2 du, 



iu)f(cj) dui. (16) 



Since hof(to) and /'(t) are Fourier transforms of each 
other, then from Parseval's theorem, 



iw/(w) 



duj 



r 2 + 7 2 <5 2 + 2 7 / 7 )/'(t)]dt (18a) 



= r 2 + 1 2 5 2 + 2 1 J[w/5(«)/j(w)]dw, (18b) 



i/'wrdt. (17) 



where we have used /(t;7) = J- 1 [/q(w)62 7 ' j ](t) and 
/'(<) = J-^iw/Co;)]^) to obtain dl8b| . Now, if the 
pulse's GDD is constant over its spectrum, /o(w) is a 
strictly real quantity, and therefore the last term on the 
RHS of (fl8b| is equal to zero, yielding (I10[) . □ 

Relation (|T0| is the reason why the ACE procedure 
can be applied for arbitrary XUV spectra (of course, 
provided that the attosecond pulse is short compared to 
the half-period of the streaking field). 
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